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Abstract: We study the generation of the baryon asymmetry in a variant of the 
standard model, where the Higgs field is stabilized by a dimension-six interaction. 
Analyzing the one-loop potential, we find a strong first order electroweak phase 
transition for Higgs masses up to at least 170 GeV. Dimension-six operators induce 
also new sources of CP violation. We compute the baryon asymmetry in the WKB 
approximation. Novel source terms in the transport equations enhance the generated 
baryon asymmetry. For a wide range of parameters the model predicts a baryon 
asymmetry close to the observed value. 
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1. Introduction 

The baryon asymmetry of the universe has recently been determined to an unprece- 
dented accuracy, 

r/B = — = (8.9±0.4) X 10"", (1.1) 
s 

by combining measurements of the cosmic microwave background [1] and large scale 
structures [2]. Explaining its origin is one of the great challenges of modern particle 
physics and cosmology. For baryogenesis Sakharov's three conditions, B violation, 
CP violation and deviation from thermal equilibrium have to be satisfied. In principle 
these conditions could be met within the standard model (SM) at the electroweak 
phase transition (EWPT) [3]. A more quantitative analysis shows however that the 
baryon asymmetry cannot be explained within the SM because there is not enough 
CP violation [4] and the phase transition turns into a smooth crossover for Higgs 
masses rriH > 80 GeV [5]. In fact, electroweak baryogenesis requires an even stronger 
criterion to be satisfied: The Higgs vacuum expectation value (vev) at the critical 
temperature, Vc = {(j){Tc)), must be larger than about Tc in order to avoid baryon 
number washout after the phase transition. 

Motivated by the possibility that electroweak baryogenesis can be tested at future 
colliders, there were many proposals to realize this mechanism in extended settings 
[6]. Some recent attempts can be found in ref. [7]. In supersymmetric models a 
strong EWPT can be induced by a light top squark [8]. Supersymmetry breaking 
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also provides new sources of CP violation. However, by now this scenario is quite 
constrained due to the negative Higgs searches. In the SM a lower bound of rriH >114 
GeV was established [10]. A strong first order phase transition could also be driven 
by cubic interactions of a singlet Higgs field [9]. 

Recently, an alternative idea caught attention: non-renormalizable operators 
could have an impact on the EWPT. These operators parametrize the effects of 
new physics beyond some cut-off scale M. In order to be relevant at weak scale 
temperatures we have to assume that M < 1 TeV. This new dynamics could be an 
ordinary quantum field theory, e.g. an extended Higgs sector. It might as well be 
something more fundamental, like strongly coupled gravity if the hierarchy problem 
is solved by the presence of extra dimensions. 

If the Higgs potential is stabilized by a 4>^ interaction, a strong first order phase 
transition can occur for Higgs masses well above 100 GeV [11-13]. A first order 
transition is triggered by a barrier in the Higgs potential. It can be provided by the 
one-loop thermal corrections of the weak gauge bosons. In the model under consid- 
eration, a barrier can also be generated from a negative 0^ term, which no longer 
destabilizes the Higgs potential. The latter possibility turns out to be dominant in 
a large part of the parameter space. Non-renormalizable interactions also allow for 
new sources of the CP violation to fuel baryogenesis [14, 15]. 

In the following we will investigate the strengh of the EWPT in the SM with low 
cut-off, taking into account the one-loop corrections to the potential. At one-loop 
the phase transition is somewhat weaker than found in the analysis of ref. [12], where 
only the thermal mass part of the one-loop correction was taken into account. Still 
we find a strong first order EWPT for Higgs masses up to at least 170 GeV, if we 
require M > 500 GeV. We will study the properties of the bubble profile, finding in 
particular that the wall thickness varies in a wide range from about 2 to 16 times 
1/Tc. We will discuss dimension-6 interactions between the Higgs field and the top 
quark which provide the necessary CP violation to fuel baryogenesis. In the WKB 
approximation these operators induce CP violating terms in the top quark dispersion 
relation which vary along the bubble wall and enter the transport equations as source 
terms. We will discuss novel source terms in the transport equations which enhance 
the generated baryon asymmetry. We find that the model can explain the observed 
baryon asymmetry for natural values of the parameters. 

2. The strength of the phase transition 

The dynamics of the EWPT is determined by the effective Higgs potential. As 
proposed in ref. [11], we add a non-renormalizable 0^ operator to the SM potential, 
so that 
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where 0^ = 2$"^$ with the SM Higgs doublet <l>. 

At finite temperature we add a thermal mass term to the potential. Because 
of the positive definite 0^-term, the quartic coupling A is allowed to take negative 
values. In the high temperature expansion of the one-loop thermal potential we get 
the thermal Higgs mass term 

^ ^\^+^9l + ^9i' + lyt)T'<l>', (2.2) 



2 V2 IQ IQ 4' 

where yt is the top Yukawa coupling and g2 and gi are the SU{2)l and U{1)y gauge 
couplings. We also include the one-loop contributions due to the transverse gauge 
bosons 

Tcf)^ (2.3) 



167r 

and the top quark 

to the effective potential, where cp ~ 13.94 [16]. We choose Q = mtop = 178 GeV. 
Another choice of Q would only change the value of the self-coupling A. Moreover 
we add the leading one-loop and two-loop corrections due to the 0^ interaction 

^ (20^r' + 0'r^). (2.5) 



8M2 

Altogether we end up with the high temperature effective potential 

T) = \ [-^.^ + Qa + ^^gl + 1,,^ + \y^) ) 0^ 



16.^^^ +4^^+ 64^^*^^^ ^^c.TV 

+ ^(06 + 20^r2 + 02r^). (2.6) 



With the two conditions 

5Kff(0,o) 



where 



<j>=v 



a^Kfr(0,O) 2 .0 7^ 



Kff(0,O)=y(0)-gA,y,Vln(^) (2.8) 

is the zero-temperature potential including the one-loop correction from the top- 
quark, we can express the two parameters ^ and A by the physical quantities mu 
and V — 246 GeV. In the following we take mn and M as the free parameters of 
the model. The SM bound on the Higgs mass apphes to our model, so we require 
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Figure 1: Contours of constant ^ = Vc/T^ in the M-m/^-plane. M and mn are given in 
units of GeV. Below the lowest line the zero-temperature minimum at 7^ is no longer 
the global one. Below the metastability line the probability for thermal tunneling gets too 
small. 

rriH > 114 GeV. We need M < 1 TeV in order for the dimension- six operator to be of 
relevance for the phase transition. If M becomes too small, for a fixed value of mn, 
the electroweak minimum ceases to be the global minimum of the zero-temperature 
potential. As shown in fig. ^ this happens around M < 500 GeV, and we exclude 
these values from the parameter space. 

During a first order phase transition there exist two energetically degenerate 
phases separated by an energy barrier at the critical temperature T^. To obtain 
and the non-zero value of the vacuum expectation value Vc the two conditions 



have to be fulfilled. The critical temperature in case of the EWPT is around 100 
GeV. At some particular temperature below Tc, say T„ (nucleation temperature), 
the broken phase bubbles nucleate, expand and percolate. The Higgs field changes 
rapidly as the bubble wall passes through space. Baryogenesis has to take place 
outside the bubble while within the bubble the sphaleron induced {B + L)-violating 
reactions must be strongly suppressed. Otherwise the generated baryon asymmetry 
would be washed out after the phase transition. The sphaleron rate is practically 
switched off if the "washout criterion" [17] 








and 




(2.9) 




> 1.1 



(2.10) 
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is satisfied. This is the condition for a first order transition to be strong. As was 
discussed in ref. [12], the sphaleron energy and therefore eq. ( p.lO|) are practically 
not affected by the presence of the 0^ term. 

In fig. |l| we show the strength of the phase transition as a function of the model 
parameters. As expected the EWPT becomes weaker for increasing Higgs masses. 
For the smallest allowed Higgs mass we need M < 825 GeV to satisfy the washout 
criterion. In contrast to the SM we find a strongly first order phase transition, even 
for Higgs masses above 150 GeV. A large part of the parameter space meets the 
requirements of electroweak baryogenesis. As M approaches the region of wrong 
zero-temperature minimum, the critical temperature becomes smaller and ^ larger. 
For ^ 3 the high temperature approximation breaks down for the top quark. 

What Higgs masses are compatible with the washout criterion depends on how 
small M is allowed to be. There is no particular bound on the (p^ operator [18]. 
However, in fig. |1| we take M > 400 GeV in order to make an expansion in powers 
of v/M reasonable. In an effective field theory all operators which are allowed by 
the symmetries are expected to be present. In particular, we expect dimension-6 
operators involving gauge fields, such as (1/M^)($"I'D^$)^. These operators have to 
be suppressed by a higher scale of about 10 TeV in order to be in agreement with 
the electroweak precision data [12]. Thus a tuning of couplings on the order of (10 
TeV/M)^ is required, and has to be explained by the UV completion of the model. 

At the one-loop level the phase transition is somewhat weaker compared to the 
analysis of ref. [12]. There only the thermal masses ( p.2| ) were included in the com- 
putation. For instance, taking rriH = 115 GeV, we find M = 825 GeV to arrive at a 
strong enough EWPT, i.e. ^ = 1.1. Including only the thermal mass corrections, one 
arrives at ^ = 1.43, and the cut-off scale can be increased to about 870 GeV until 
the phase transition becomes too weak [12]. 

How important are the different one-loop contributions? For ^ > 1 the cubic 
term ( |2.3|) is still relevant: Leaving it out considerably weakens the phase transition 
from ^ = 1 to ^ = 0.56, for mn = 115 GeV. Omitting also the log-term ( |2.4|) makes 
the phase transition stronger again, ^ = 0.81. In addition, getting rid of the one-loop 
term of eq. ( p.5| ) increases C, to 1.27. Thus the one-loop contributions in ( |2.3| ) - ( |2.5| ) 
partially cancel each other and therefore our results agree reasonably well with those 
of ref. [12]. For larger Higgs masses and stronger phase transition the picture is 
qualitatively similar, however, the cubic term becomes less important. 

The two-loop 0^ term of eq. (|2.5|) practically does not change the result, demon- 
strating that the dimension-6 operator does not spoil the loop expansion. We have 
also checked that adding a dimension-8 term (1/M'^)($"''$)^ only affects ^ at the 
order of v^/M^. 

The one-loop effective potential was also discussed in ref. [13]. However, the 
authors impose an erroneous lower bound on the cut-off scale, requiring a positive 
mass squared for the Goldstone boson. As a result, they obtain a lower bound on 
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the Higgs mass from eq. ( 2.10 ) which is much too small. 



3. The bubble properties 

In this section we discuss some bubble properties which will enter the computation 
of the baryon asymmetry, in particular the thickness L„, and the velocity, of the 
wall. As already mentioned, the two minima of become of the same depth at T^, 
but tunneling with the formation of bubbles of the broken phase will start somewhat 
later at a temperature T„. The probability for thermal tunneling is exponentially 
suppressed by the energy of the critical bubble, S'3. The phase transition starts if the 
nucleation probability per horizon volume becomes of order unity, which translates 
to ^3(T„)/T„ - 130 - 140 [19]. 

For ^ = 1 the amount of supercooling, i.e. the difference between the critical and 
nucleation temperatures, is small. For mu = 115 GeV we find = 107.34 GeV and 
Tc — Tn = 0.45 GeV. The system is well described by the thin wall approximation. 
For smaller values of M and stronger phase transition, supercooling becomes more 
and more important. The thin wall approximation is no longer reliable. At some 
critical value the phase transition does no longer proceed at all. The universe remains 
stuck in the symmetric vacuum. This regime is indicated by the "metastability" line 
in fig. g. 

Once a critical bubble is nucleated it will expand. The expansion is accelerated 
by the internal pressure and slowed down by plasma friction. Finally, a stationary 
situation will be reached, where the different forces are balanced, and the wall prop- 
agates with constant velocity, v^. In order to estimate the thickness of the bubble 
wall, we ignore friction for a moment and solve the field equation at the critical 
temperature with the effective potential of eq. (p.6|). 



The boundary conditions read (f){z — > —00) = Vc and (f){z — > 00) = 0. The bubble 
profile can approximately be described by a kink, 

cPiz) = 1(1- tanh-^^ (3.2) 

with = \/vfJ{8Vb), where H is the height of the potential barrier. This relation 
would be exact for a 0^ potential and we found that it is also a good approximation 
in our case. In fig. ^ the wall thickness is shown in dependence of the Higgs mass 
rriH and M. As we decrease M at fixed rriH, the wall thickness in units of 1/Tc 
becomes smaller. The same happens if we decrease ttlh at fixed M. The main effect 
comes from the decrease of Tc in these cases. Notice that L^^Tc varies in a wide range 
between about 2 and 16. 
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Figure 2: The wall thickness L„ as a function of the Higgs-mass rriH for several values of 
M, which are given in units of GeV. In addition lines of constant ^ are shown. 



Let us finally comment on the wall velocity. Taking into account only the friction 
related to the infrared gauge field modes [20] , 

^ 327rL„ AV^ 

we obtain a wall velocity of order unity, except for being close to one. Here AV^ 
is the potential difference at the nucleation temperature and mw the mass of the W 
boson. The order unity correction in the denominator is induced by friction of other 
particles in the plasma, in particular the top quark [21]. Since numerically \ia{mwL^) 
is only of order unity, top quark friction will slow down the wall considerably. The 
wall velocity is further reduced by latent heat of the nucleating bubbles. In general, 
the wall moves faster in the case of a stronger phase transition. 

Let us briefly discuss two representative examples. Taking ^ = 1 and mn = 115 
GeV, we obtain = 0.24 from eq. (^.31). Including the effect of reheating, AV is 
reduced and we arrive at = 0.08. If we finally switch on the order unity correction, 
we end up with ~ 0.05. The picture looks very much different for stronger phase 
transitions. Going to ^ = 1.5, eq. ( p.3|) already leads to a wall velocity of order unity. 
Including again top quark friction and reheating, we obtain ~ 0.5. For larger 
Higgs masses we find a very similar behavior. These are only very rough estimates, 
since eq. ( |3.3| ) breaks down for large values of f„. Given these uncertainties we will 
treat clS db free parameter in our computation of the baryon asymmetry. 
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4. CP violation 



Non-renormalizable interactions provide new sources of CP violation. In the absence 
of gauge singlets the leading operators are of mass dimension six. In ref. [14] a 
operator was discussed. We will focus on the operators 



^J$t$)$„c^^. (4.1) 

which have been proposed to drive baryon number generation in ref. [15]. There are 
analogous terms for the down-type quarks and leptons. The fermion masses become 

V V 

=l/..^ + ^^7|^^.. (4.2) 

where Uij are the ordinary Yukawa couplings. In unitary gauge, the effective Yukawa 
couplings to the physical Higgs boson are given by 

3^2 

Thus, there is a mismatch of order {xv'^/M'^) between the fermion masses and the 
effective Yukawa couplings. In general, the couplings Xij contain complex phases, and 
they are of unknown flavor structure. While for the top quark Xgg may be of order 
unity, the couplings of the lighter fermions should not exceed 0{M'^mf/v^) to avoid 
fine tuning of the fermion masses ruf. For instance, the corresponding coupling for 
the electron should at most be of order 10~^ x (M/TeV)^, which is a small number. 
Having this in mind, and lacking a theory of flavor mixing, we will therefore assume 
that the Xij have a similar flavor structure as the corresponding Yukawa couplings, 
i.e. Xij ~ Uij, up to order unity coefficients. This structure could be motivated by 
a Froggatt-Nielsen [22] type mechanism, where the operators ( [4.1| ) and the ordinary 
Yukawa couplings would have the same quantum numbers. 



Since the operators ( |4.1| ) do not have to be strictly aligned with the ordinary 
Yukawa couplings, they will induce flavor violation and CP violation. Via tree-level 
Higgs exchange they generate four-fermion interactions, which, for instance, affect 
K — K mixing. For example, the operator {1/Q'^){d'^sds'^) is generated at a level of 

^ -42{4ir- (4.4) 



Q2 M^m]j 

For ~ sin6'c/t', where denotes the strange quark mass and dc the Cabibbo 
angle, we obtain roughly Q ~ 5 ■ lO^GeVx (M/TeV)^. At M > 30 TeV a 2-loop 
contribution of order {1 / {lQTt'^Y)xi2{x2i)* / M'^ becomes dominant.^ Experimentally, 
these operators are constrained to be suppressed by Q > 10"^ GeV, especially in the 



^This diagram is quadratically divergent which we cut off at the scale M . 



8 



presence of CP violation. This constraint is compatible with (|4.4|) for M > 500 
GeV. Later on we will only be interested in the coupling of the top, X33 ~ 1. If we 
make the unnatural assumption that only this coupling is present, a ddss operator 
is generated at 2-loops with Q ~ 4 ■ lO^GeVx (M/TeV)^. In the B system the 
experimental bounds are even more easily satisfied. 

CP violating couplings induce electric dipole moments (EDM's). The EDM of 
the neutron, for instance, is experimentally constrained by dn/e < 6.3- 10~^^ cm [23]. 
The individual EDM's of up and down quarks should therefore be not much larger. 
The up quark receives the larger contribution. At one-loop we find 



e 167r2 V ^ / 



(4.5) 



which might be close to the experimental bound. In the second step we assumed a 
maximal phase and |x"3| ~ ~ Vub- If only the coupling X33 is present, an EDM 
of roughly du/e ~ 5 ■ 10"^'^ cm x Im(x33) x (TeV/M)^ is induced at the 2-loop level 
[15]. Thus a coupling X33 of order unity with a large phase can be present without 
inducing too large flavor changing neutral currents or EDM's, even for M ~ 300 
GeV. Other couplings also are allowed, as long as some Yukawa-like hierarchy is 
respected. However, experimental signals of non-standard flavor physics could be 
detected in the near future. 

At low energies the non-renormalizable operators and the Yukawa couplings melt 
into the couplings ([4.3|) . However, during the EWPT the two terms in (|4.3| ) vary 
differently along the bubble wall. As a result, the fermion masses acquire position 
dependent phases which cannot be rotated away. For the phase of the top quark 
mass we obtain 

tan^,(z)^sin^,^^, (4.6) 

where we defined xte^'^* = X33 and ignored the real part of X33. In two Higgs doublet 
models such a phase may arise from spontaneous CP violation. In supersymmetric 
models position dependent phases are induced by flavor mixing, e.g. for the charginos. 
In the next section we discuss how the phase ([4.6|) drives the generation of a baryon 
asymmetry as a bubble wall moves through the plasma. 



5. Transport equations 

The CP violating interactions of particles in the plasma with the bubble wall create an 
excess of left-handed quarks over the corresponding anti-quarks. ^ In the symmetric 
phase the left-handed quark density biases the sphaleron transitions to generate a 
net baryon asymmetry. 

^Quarks will turn out to be more important than leptons because of the large top mass. 
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In the WKB approach the CP violating interaction of a fermion with the wall 
leads to different dispersion relations for particles and anti-particles [24], depending 
on their complex masses. To make this method applicable, it is required that the 
typical de Broglie wavelength of particles in the plasma is small compared to the 
width of the wall, i.e. L„T ^ 1. Otherwise an expansion in derivatives of the 
background Higgs field cannot be justified. According to the results of section 3 this 
is a good approximation in a large fraction of our parameter space. It is violated only 
in the cases of a very strong phase transition, ^ > 3. From the dispersion relations 
we compute a force term which enters the transport equations that describe the 
evolution of the plasma. An alternative approach was followed in ref. [26]. 

Let us consider a single Dirac fermion, such as the top quark, with a space-time 
dependent mass ReA4{z) + i'y^lmA4{z), where Ai{z) = m{z)e^^^^^ and z = z — v^t 
denotes the relative coordinate perpendicular to the wall. For particles and anti- 
particles the dispersion relations to first order in derivatives read [25] 



2 

m 



E± = Eo±AE= v/p2 + m2±sign(p,)e'— — -, (5.1) 

2(p^ -|- m^) 

where = is the squared kinetic momentum and 6' = dO/dz. is the energy 
of left-handed particles, corresponds to the right-handed states, and for the anti- 
particles the other way round. ^ In a more rigorous treatment similar dispersion 
relations were derived for spin states in the Schwinger-Keldysh formalism [27, 28] . 

The evolution of the particle distributions /j(t,x, p) we describe by classical 
Boltzmann equations. The dispersion relations ( |5.1| ) induce force terms, which are 
different for particles and anti-particles. To make the system of equations tractable, 
we use a fluid-type ansatz in the rest frame of the plasma [24] 

/.(t,x,p) = ^^^^^_J^_^^^^^ , (5.2) 

where f , and /Xj denote velocity perturbations and chemical potentials for each fluid. 
The velocity perturbations are introduced to model the movement of particles in 
response to the force. 

For a shorter notation let us first introduce some symbols K, which represent 
momentum averages normalized relative to the massless Fermi-Dirac case, 

^ Jd^pXjp) 
^^^-/d3p/;(m = 0)' ^^-^^ 



^Later on the relevant particles will be relativistic, so that we can approximate helicity by 
chirality. 
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where /4 = -/^e^^Vle^'^" ± 1)^- We define 




Kr, = ( J'%,, mm)), (5.4) 



which appear in the transport equations discussed in the following. In the case of 
a massless fermion we obtain Ki = 1.096T, K2 = A.QOQT'^, K3 = 0.211/T, = 
0.105/T, = 0.105/T, Kq = -0.038/T3 and K7 = O.IOS/T. For m > T the 
averages experience an exponential Boltzmann suppression. 

We look for solutions of the Boltzmann equation which are stationary, i.e. which 
only depend on the relative coordinate z. We expand the Boltzmann equation in 
derivatives of the fermion mass. At first order in derivatives there is no difference 
between particles and anti-particles. Weighting the Boltzmann equation with 1 and 
Pz, we obtain after momentum averaging 

- Ki^iv[i - {Ci) = Ks^iV^^mi)' (5.5) 
- Ki,ifi[ .^ + K2,iV„vl^ - {pzCi) = 0. (5.6) 

Here /Xj^i and fj 1 indicate the perturbations to first order in derivatives. A prime 
denotes again a derivative with respect to z. The statistical factor Kj = {f'^irrii)) is 
1 (2) for massless fermions (bosons) and becomes exponentially small for m ^ T. 
The force term on the right-hand side is induced by the change in the particle mass 
along the wall. Introducing inelastic rates, r™'^^ and elastic rates, , for a process 
p, the collision terms take the form [29] 

{Q = Yl E ' (P^^^) = E = ^^p'rf . (5.7) 



In {pzCi) we neglected inelastic processes. The leading order eqs. (|5.5|), ( |5.6| ) contain 
the friction terms which enter the computation of the wall velocity [21]. 

To second order in derivatives, we have to distinguish between particles and 
anti-particles. The perturbations contain CP odd and even components, 

fJ-i = fJ'i,! + f^i,2o + fJ'i,2e, Vi = Vi^i + Vi^2o + Vi,2e- (5.8) 
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In the following only the odd second order perturbations will enter, so we can drop 
the subscript "o" to simplify the notation. Subtracting the equations of particles 
and anti-particles, we obtain 

l^iVwlJ''L2 - - {Q = -^6,i^i"^iyLti,l (5.9) 

--^i,iK,2 + K2,iV^v[ 2 - {PzCi) = Ki^iV^mje'l + K^^iV^{'ml)'e[ - K'j^imje[v[ ,^.{h.lQ) 

Note that the CP violating source terms are proportional to derivatives of ^^j. A 
constant phase does not contribute. The source terms proportional to the first order 
perturbations have not been investigated so far in a realistic context. (See ref. [28] 
for a discussion in the context of Schwinger-Keldysh formalism.) To study their 
relevance for the generation of the observed baryon asymmetry will be a main issue 
in the next section. 

We may use eq. (|5.9|) to solve for fj 2- Neglecting derivatives of the thermal 
averages, which are higher order in derivatives, we end up with diffusion equations 
for the chemical potentials 

-K, A(l - AivDn'l^ - K^v^f^'i^2 - AAww 

P 3 

+ E E /^^-.^ - A At^w E(rr')' E ^^-^ = (^-n) 

P 3 P 3 

where 

A = (5.12) 
In the massless limit we have A ^ 3.83. The diffusion constants are given by [29] 

'^^A = (5.13) 

To leading order in the wall velocity, neglecting derivatives of the inelastic rates and 
ratios of inelastic to elastic rates, the left-hand side of eq. ( ^.11| ) reproduces the result 
of ref. [29]. This corresponds to dropping the terms proportional to Ai. In the next 
section we will examine to what extent this simplification is justified. Since Ai is not 
a small number, the corrections will turn out to be important in certain regimes. For 
the source terms we obtain 

= {K,,m^e': + K,Amm' 

(i - ^At^w^) K^:/.: J - ^^^^ K^:^a)' . (s.m) 

In the second line new source terms related to the first order perturbations are 
showing up. As expected, no source terms are left in the case of vanishing wall 
velocity. 
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Let us apply these general results to the SM with a low cut-off. In a first 
step we compute the asymmetry in the left-handed quark density. At this stage 
we neglect the weak sphalerons, i.e. baryon and lepton number are conserved. The 
most important particle species are the left- and right-handed top quarks and the 
Higgs bosons. Leptons are only produced by small Yukawa couplings and therefore 
not taken into account. It turns out that also the Higgs bosons have only a minor 
impact on the generated baryon asymmetry. They change the final result only at the 
percent level, so we can ignore them. In a second step, the weak sphalerons convert 
the left-handed quark number into a baryon asymmetry. 

We take into account the top Yukawa interaction, Ty, the strong sphalerons, Tgs 
and the top helicity flips, caused by the interactions with the bubble wall. The 
latter are only present in the broken phase. The gauge interactions are assumed to 
be in equilibrium. The transport equations become 

-6r„ (/ig3,2 + /it,2) - Q^ss [(2 + 9Kt + 9Kb) /ig3,2 + (1 - 9/tt) /it,2] 

= -3ir6X^K,i (5-15) 



2 

> \ 

"■'73,2 



= K,,tV„mX + K5,tv^{m^M - K7,tmXv',,,i (5.16) 

-STss [(2 + 9^4 + 9K,b) f^q3,2 + (1 - 9Kt) f^t,2] 

= -3K,,m'M^, (5.17) 



-^l,tl^'t 2 + K2,tVwV't 2 7r'"t,2 

KtDg 

= K^,tV^mX + K^,tvAm]ye[ - K^^tmle[v[^^. (5.18) 

The top quark phase, 6t, is given by eq. (|46|). For the chemical potentials we take 
= /^(^s) and /igg = (/i(M3) + ii{d^))/2. The index t {b) refers to the top and 
bottom quark, respectively. We have omitted the tiny source of the bottom quark 
which is suppressed by {mb/mtY. We used baryon number conservation to express 
the strong sphaleron interaction in terms of Hq^^2 and jj,t,2 [30]. Replacing the second 
order source terms by the first order ones, the same system of transport equation 
holds for the first order perturbations. 

Using again baryon number conservation, the chemical potential of left-handed 
quarks, /xs^ = /ig^ + jJiq^ + A^gg, is obtained as 

IIBl = (1 + 2/«t + 2K,,)/ig3 - (5.19) 
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The baryon asymmetry is then given by [25] 



oo 



ub 405r^^ , 

where f^^ is the weak sphaleron rate and u = 4:5Tws/{4:V^T^). The effective number 
of degrees of freedom in the plasma is = 106.75. In eq. ( |5.20|) the weak sphaleron 



rate has been suddenly switched off in the broken phase, z < 0. The exponential 
factor in the integrand accounts for the relaxation of the baryon number if the wall 
moves very slowly. 



6. Numerical Results 



In this section we present our evaluations of the transport equations ( p.l5[ ) - ( ]5.18| ). 



We will discuss under what conditions the terms proportional to Ai may be neglected 
in eq. (|5.11| ) and investigate what is the impact of the new source terms in eqs. (|5.9| ) 
and ( ^.10|) . Finally, we will demonstrate that the SM with a low cut-off can explain 
the observed baryon asymmetry for natural values of the parameters. 

In our numerical computations we use the following values for the weak sphaleron 
rate [31], the strong sphaleron rate [32], the top Yukawa rate [30], the top helicity 
flip rate and the quark diffusion constant [30] 

f = 1.0 ■ 10-^T^ f = 4.9 ■ 10-^T^ 

r, = 4.2 . lo-^T, = 

Dq = ^. (6.1) 

Note that in F^s = F^^T^ the value as = 0.086 from the dimensionally reduced 
theory has been used [32]. Changing the rates F^ and F^ has only a small effect 
on the baryon asymmetry, as would have the inclusion of the Higgs field chemical 
potential in the transport equations. Doubling the value of Dq enhances the baryon 
asymmetry by 20-30% because of more efficient diffusion. Enhancing Tss reduces the 
baryon asymmetry since the strong sphalerons drive (Ibl to zero if the quarks are 
taken massless [33]. The baryon asymmetry also depends on whether we take the 
top quark to be massive or massless in the thermal averages. If we switch on the top 
mass, the baryon asymmetry becomes smaller since the thermal averages go down. 
In the evaluations we use the half of the broken phase to compute the averages. 



Let us first discuss under which conditions the Ai corrections in eq. ( |5.11|) become 
important. At this stage we do not yet relate the bubble wall parameters to the model 
introduced in section 2. In fig. ^ we display the baryon asymmetry computed with 
the simplified equations (dashed lines) compared to the one obtained from the full 
equations ( p.l5| ) - ( |5.18[ ) (solid lines) as a function of the wall velocity. We take 
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Figure 3: Comparison between eq. ( |5.11| ), where the terms proportional to Ai have been 
neglected (dashed), and eqs. ( 5.15| ) - (|5.18| ) without the first order perturbations (solid) for 
different values of L„. The other parameters are taken as ^ = 1.5 and M = 6T. 



^ = 1.5 and M/T = 6 for three different values of L^. The CP violating phase in 
the dimension-6 operator ( |4.1j ) we take as maximal, i.e. simpt = 1, and we choose 
Xt = 1. In any case the baryon asymmetry is proportional to Xtsimff Since we 
want to compare the left hand sides of the transport equations, we included only the 
source term of the first line of eq. ( |5.14| ). The simplified equations give a reasonable 
description for < 0.1. For large values of there are sizable deviations, especially 
for thinner bubble walls. This behavior is expected since the Ai corrections come 
with additional powers of the wall velocity. In the MSSM, where fw ~ 0.05 — 0.1 [34], 
the simplified equations are applicable. In the following we will use the full equations 
( |5.15|) - ( |5.18|) to compute the baryon asymmetry. Fig. ^ also demonstrates that the 
baryon asymmetry increases for thinner bubble walls. This behavior is expected 
since the source terms involve derivatives of the background Higgs field. 

In fig. ^ we compare the contributions to the baryon asymmetry due to the 
different source terms on the right hand side of eqs. ( |5.9| ) and ( |5.10|) , using the 
parameters of fig. ^ and L„ = 8/T. The new source terms proportional to the 
first order perturbations Hi^i and Vi^i are non-negligible. They enhance the baryon 
asymmetry, especially for small values of v^,. For large wall velocities they do no 
longer matter. The total baryon asymmetry depends only mildly on v^, which is 
quite positive, given our poor understanding of this parameter. For other wall widths 
the picture is similar. 

As shown in fig. ^ the baryon asymmetry increases rapidly for larger values of 
^. We fixed f„ = 0.01 and 0.3 and again L„ = 8/T. For large ^ the top quark mass 
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Figure 4: The solid line represents r/s as a function of the wall velocity for = 8/T, 
M = 6T and ^ = 1.5. The dashed line would be the asymmetry without the source terms 
containing the first order perturbations n'l and v'l and the dotted one is the contribution 
due to these terms only. 



becomes larger in the broken phase, and the source terms involve powers of mf. Also 
the CP violating phase in the top quark mass from eq. ([4.6|) becomes stronger. The 



source terms related to the first order perturbations have an additional power of mf. 
Therefore they grow even faster and dominate for large ^. This behavior holds also 
for other values of L„. 

Let us finally relate the bubble wall parameters to our SM with low cut-off. In 
the fig. P we present the baryon asymmetry in the model for rriH = 115 and 150 GeV 
as a function of the cut-off scale M. For every value of M we compute the strength of 
the phase transition and the bubble width. We consider one small = 0.01 and one 
large wall velocity fw = 0.3. We take again a maximal CP violating phase sin (ft = 1 
and Xt = 1. As expected the baryon asymmetry grows rapidly as we lower M. At the 
very lowest values of M the wall thickness becomes of order 1/T (see fig. ^ so that 
our WKB approach ceases to be reliable. Moreover, the bubble walls may become 
relativistic in this regime, and diffusion of charges into the symmetric phase may no 
longer be efficient. One can see that rjE depends only mildly on the wall velocity. 

Nevertheless, independent of the Higgs mass we have chosen, we can generate the 
observed baryon asymmetry (|1 . 1|) without amplifying the CP violating dimension-6 
operator (4J^). This requires the phase transition to be sufficiently strong, i.e. ^ > 1.5. 



At this strength of the phase transition our computation of the baryon asymmetry 
is still under control. For smaller values of ^ we have to take Xt > 1. For instance, in 
the case of ttih = 115 GeV and ,^ = 1.1 (M = 825 GeV), we could use Xt ~ 40 which 
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Figure 5: The baryon asymmetry as a function of ^ for two wall velocities. The different 
line types have the same meaning as in fig. ^, and again L„ = 8/T and M = 6T. On the 
right edge of the figure the upper curves are for t;^ = 0.01 and the lower ones for = 0.3. 




500 550 600 650 700 750 800 850 

M 

Figure 6: The baryon asymmetry in the SM with low cut-off for two different Higgs masses 
as a function of M (in units of GeV) for = 0.01 (solid) and = 0.3 (dashed). The 
horizontal lines indicate the errorband of the measured value. 

is barely consistant with the bound from the neutron EDM discussed in section 3. 
Thus the SM with low cut-off can account for the observed baryon asymmetry in a 
wide range of the model parameters, without being in conflict with constraints from 
flavor and CP violation. 
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At the end of this section let us briefly comment on the impact of sources from 
the bottom quark and the tau lepton. In case of the bottom the source term is 
heavily suppressed by {rrib/mt)^ ~ 10^''. The tau lepton is more relevant because of 
its larger diffusion constant of about 380/T [35]. Still its contribution is about 10^ 
times smaller than that of the top quark and can be safely neglected. 

7. Conclusions 

We have investigated the electroweak phase transition and baryogenesis in the stan- 
dard model augmented by a dimension-6 Higgs self interaction. Taking the sup- 
pression scale of this operator to be M < 1 TeV, the EWPT becomes first order, 
without introducing new degrees of freedom in the model. In addition to the Higgs 
mass only the parameter M enters the computation of the phase transition. There is 
no relevant bound on the (j)^ interaction. However, dimension-6 operators involving 
for instance gauge fields, which are also expected to be present in a general effective 
field theory, have to be tuned at the level of 10~^ in order not to spoil the electroweak 
fit [12]. 

Requiring M > 500 GeV the phase transition is strong enough to prevent baryon 
number washout for Higgs masses up to 170 GeV. In our analysis we have used the 
one-loop thermal potential. The phase transition is somewhat weaker than found in 
ref. [12], where only the thermal mass part of the one- loop correction was taken into 
account. We have checked that the dimension-6 operator does not spoil the loop 
expansion of the effective potential. We have computed the wall thickness which 
turns out to vary in a wide range from about 2/T to 16/T. As M becomes smaller 
the EWPT becomes stronger and the bubble wall thinner. For very small cut-off 
scales the symmetric phase becomes metastable. 

A dimension-6 operator involving the Higgs field and the top quark provides a 
new source of CP violation. It induces a complex phase in the top quark mass which 
varies along the phase boundary. We discussed that this operator is consistent with 
present bounds on EDM's and flavor violation for M > 300 GeV. However, it may 
leave a detectable signal in forthcoming experiments. 

As a result of the varying phase, top quarks and anti-top quarks experience a 
different force as they pass through the bubble wall. We treat the system in the 
WKB approximation, expanding in derivatives of the background Higgs fleld. Our 
approach is valid for a large fraction of the parameter space of the model. It will break 
down for very small values of the cut-off scale M, where the bubble width becomes of 
order the inverse temperature. The CP violating force enters the transport equations 
which describe the hot plasma. Carefully expanding in derivatives of the wall proflle, 
we flnd novel source terms which enhance the generated baryon asymmetry. They 
are especially relevant for slow bubble walls and dominate over the known source 



18 



terms for large values of the particle mass. Because of these properties they should 
play a prominent role in the MSSM, where the top quark is replaced by the charginos. 

In the model considered, the observed baryon asymmetry can be explained for 
natural values of the parameters. The phase transition should be somewhat stronger 
than required by the washout criterion. If the EWPT is not that strong, the coeffi- 
cient of the CP violating dimension-6 operator has to be taken larger than one, which 
is compatible with experiments. It would also be interesting to study the impact of 
other CP violating operators, such as the one discussed in ref. [14], which has been 
ignored in our study. 

With a low cut-off the model is expected to lead to non-standard signals in flavor 
physics, such as EDM's and flavor changing neutral currents, which can be tested in 
future experiments. The LHC will be able to directly test the physics at the cut-off 
scale. If the general cut-off scale is in the multi-TeV range and the 0^ interaction is 
anomalously large, the model could still be identified by its non-standard Higgs self 
couplings. However, the required precision will probably take a linear collider [12]. 

In conclusion, the standard model with low cut-off provides the missing ingre- 
dients for electroweak baryogenesis: a strong phase transition and additional CP 
violation. Moreover, its simple structure makes it an ideal laboratory to refine the 
computation of the baryon asymmetry. 

Acknowledgements 

We thank Dominik Stockinger, Steffen Weinstock and Oscar Vives for helpful dicus- 
sions. The work of D.B., L.F., and M.S. was supported by the DFG, grant FOR 
339/2-1. 

References 

[1] D.N. Spergel et al, WMAP Collaboration, Astrophys. J. Suppl. 148 (2003) 175 
[astro-ph/0302209]. 

[2] M. Tegmark et al, SDSS Collaboration, Phys. Rev. D69 (2004) 103501 
[astro-ph/0310723]. 

[3] V. A. Kuzmin, V. A. Rubakov and M. E. Shaposhnikov, Phys. Lett. B155 (1985) 36. 

[4] M.B. Gavela, P. Hernandez, J. Orloff, O. Pane, C. Quimbay, Nucl. Phys. B430 
(1994) 382 [hep-ph/9406289]; P. Huet, E. Sather, Phys. Rev. D51 (1995) 379 
[hep-ph/9404302]; T. Konstandin, T. Prokopec and M.G. Schmidt, Nucl. Phys. 
B679 (2004) 246 [hep-ph/0309291]. 

[5] K. Kajantie, M. Laine, K. Rummukainen, M.E. Shaposhnikov, Phys. Rev. Lett. 77 
(1996) 2887 [hep-ph/9605288], Nucl. Phys. B493 (1997) 413 [hep-lat/9612006]; 
F. Csikor, Z. Fodor and J. Heitger, Phys. Rev. Lett. 82 (1999) 21 [hep-ph/9809291]. 



19 



[6] For a review, see: A. Riotto, M. Trodden, Ann. Rev. Nucl. Part. Sci. 49 (1999) 35 
[hep-ph/9901362]. 

[7] S. Kanemura, Y. Okada, E.Scnaha, hcp-ph/0411354; S.W. Ham, Y.S. Jcong, S.K. 
Oh, hcp-ph/0411352; S.W. Ham, S.K. Oh, D. Son, hep-ph/0411012; M. Carena, 
A. Megevand, M. Quiros, C.E.M. Wagner, hep-ph/0410352; A. Menon, 
D.E. Morrissey, C.E.M. Wagner, Phys. Rev. D70 (2004) 035005 [hep-ph/0404184]. 

[8] M. Carena, M. Quiros, C.E.M. Wagner, Phys. Lett. B380 (1996) 81-91 

[hep-ph/9603420]; D. Bodeker, P. John, M. Laine, M.G. Schmidt, Nucl. Phys. B497 
(1997) 387 [hep-ph/9612364]; B. de Carlos, J.R. Espinosa, Nucl. Phys. B503 (1997) 
24 [hep-ph/9703212]; M. Laine, K. Rummukainen, Nucl. Phys. B535 (1998) 423 
[hep-lat/9804019]. 

[9] M. Pietroni, Nucl. Phys. B402 (1993) 27 [hep-ph/9207227]; A.T. Davies, 

CD. Proggatt, R.G. Moorhouse, Phys. Lett. B372 (1996) 88 [hep-ph/9603388]; 
S.J. Huber, M.G. Schmidt, Eur. Phys. J. CIO (1999) 473 [hep-ph/9809506]; 
S.J. Huber, M.G. Schmidt, Nucl. Phys. B606 (2001) 183 [hcp-ph/0003122]; 
M. Bastero-Gil, C. Hugonic, S.F. King, D.P. Roy, S. Vempati, Phys. Lett. B489 
(2000) 359 [hep-ph/0006198]. 

[10] LEP Collaborations, LEP Electroweak Working Group, SLD Electroweak Group and 
SLD Heavy Flavor Group, hep-ex/0312023. 

[11] X. Zhang, Phys. Rev. D47 (1993) 3065 [hep-ph/9301277]. 

[12] C. Grojean, G. Servant and J.D. Wells, hep-ph/0407019. 

[13] S.W. Ham and S.K. Oh, Phys. Rev. D70 (2004) 093007 [hep-ph/0408324]. 

[14] M. Dine, P. Huet, R. Singleton and L. Susskind Phys. Lett. B257 (1991) 351. 

[15] X. Zhang, S.K. Lee, K. Whisnant and B.L. Young, Phys. Rev. D50 (1994) 7042 
[hep-ph/9407259]. 

[16] L. Dolan, R. Jackiw, Phys. Rev. D9 (1974) 3320. 

[17] G. D. Moore, Phys. Rev. D59 (1999) 014503 [hep-ph/9805264]. 

[18] V. Barger, T. Han, P. Langacker, B. McElrath, P. Zerwas, Phys. Rev. D67 (2003) 
115001 [hep-ph/0301097]. 

[19] G.W. Anderson, L.J. Hah, Phys. Rev. D45 (1992) 2685; for a non-perturbative 

treatment of bubble nucleation, see G. D. Moore and K. Rummukainen, Phys. Rev. 
D63 (2001) 045002 [hep-ph/0009132]. 

[20] G.D. Moore, JHEP 0003 (2000) 006 [hep-ph/0001274]. 

[21] G.D. Moore, T. Prokopec, Phys. Rev. D52 (1995) 7182 [hep-ph/9506475]. 



20 



[22] C. D. Proggatt and H. B. Nielsen, Nucl. Phys. B147 (1979) 277. 
[23] P.G. Harris et al., Phys. Rev. Lett. 82 (1999) 904. 

[24] M. Joyce, T. Prokopec and N. Turok, Phys. Rev. Lett. 75 (1995) 1695, Erratum-ibid. 
75 (1995) 3375 [hep-ph/9408339]; Phys. Rev. D53 (1996) 2958 [hep-ph/9410282]. 

[25] J.M. Cline, K. Kainulainen, Phys. Rev. Lett. 85 (2000) 5519 [hep-ph/0002272]; 
J.M. Cline, M. Joyce and Kimmo Kainulainen, JHEP 0007 (2000) 018 
[hep-ph/0006119]. Erratum hep-ph/0 110031. 

[26] M. Carcna,, J.M. Moreno, M. Quiros and C.E.M. Wagner, Nucl. Phys. B599 (2001) 
158 [hep-ph/0011055]; M. Carena, M. Quiro, M. Seco and C.E.M. Wagner, Nucl. 
Phys. B650 (2003) 24 [hep-ph/0208043]. 

[27] K. Kainulainen, T. Prokopec, M.G. Schmidt and S. Weinstock, JHEP 0106 (2001) 
031 [hep-ph/0105295], Phys. Rev. D66 (2002) 043502 [hep-ph/0202177]. 

[28] T. Prokopec, M.G. Schmidt and S. Weinstock, Annals Phys. 314 (2004) 208 
[hep-ph/0312110]. Annals Phys. 314 (2004) 267 [hep-ph/0406140]. 

[29] J.M. Chne, M. Joyce and K. Kainulainen, Phys. Lett. B417 (1998) 79, Erratum-ibid. 
B448 (1999) 321 [hep-ph/9708393]. 

[30] P. Huet, A.E. Nelson, Phys. Rev. D53 (1996) 4578 [hep-ph/9506477]. 

[31] G.D. Moore, Phys. Rev. D62 (2000) 085011 [hep-ph/0001216]; hep-ph/0009161; 
G. D. Moore, Nucl. Phys. B568 (2000) 367 [hep-ph/9810313]; P. Arnold and 
L. G. Yaffe, Phys. Rev. D62 (2000) 125014 [hep-ph/9912306]; D. Bodeker, Phys. 
Lett. B426 (1998) 351 [hep-ph/9801430]. 

[32] G.D. Moore, Phys. Lett. B412 (1997) 359 [hep-ph/9705248]. 

[33] G.F. Giudice, M.E. Shaposhnikov, Phys. Lett. B326 (1994) 118 [hep-ph/9311367]. 

[34] P. John, M.G. Schmidt, Nucl. Phys. B598 (2001) 291, Erratum ibid. B648 (2003) 
449 [hep-ph/0002050]. 

[35] M. Joyce, T. Prokopec, N. Turok, Phys. Rev. D53 (1996) 2930 [hep-ph/9410281]. 



21 



